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Abstract 



We report density-functional studies of several small molecules 
(H2, N%, CO, H2O, and C-H4) within the Krieger-Li-Iafrate (KLI) approx- 
imation to the exact Kohn-Sham local exchange potential, using a three- 
dimensional real-space finite-difference pseudopotential method. It is found 
that exchange-only KLI leads to markedly improved eigenvalue spectra com- 
pared to those obtained within the standard local-density approximation 
(LDA), the generalized gradient approximation (GGA), and the Hartree-Fock 
(HF) method. For structural properties, exchange-only KLI results are close 
to the corresponding HF values. We find that the addition of LDA or GGA 
correlation energy functionals to the KLI exact exchange energy functional 
does not lead to systematic improvements. 

PACS numbers: 02.70Bf, 31.10.+Z, 31.15.Ew, 31.15.Fx, 31.50,+w, 31.90+s 
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I. INTRODUCTION 



Density functional theory (DFT) within the Kohn-Sham (KS) single-particle formalism 
has now been established as a standard tool for the study of inhomogeneous quantum many- 
electron systems. [l|-|3| The local density approximation (LDA) and generalized gradient 
approximation (GGA) are successful in describing the ground state properties of various 
systems. On the other hand, their eigenvalue spectra are known to be far from the reality, 
and it has been assumed that KS orbital eigenvalues are auxiliary quantities without any 
physical meaning except the highest occupied one. 0] However, it is clear that DFT is also 
a legitimate route to the study of excited state properties via perturbation theory [Q or 
time-dependent DFT. [5]|| Especially, Gorling has recently formulated a KS formalism for 
excited states based on DFT perturbation theory and demonstrated that the KS eigenvalue 
difference is a well-defined zeroth order approximation to excitation energies. This has 
been also explicitly confirmed for several atoms by accurate quantum Monte-Carlo (QMC) 
calculations. ||||] 

It appears that - at least for localized systems - the unrealistic description of eigenvalue 
spectra is mainly due to the deficiencies of the conventional approximations for the exchange- 
correlation energy functional such as the LDA or GGA rather than the KS scheme itself. 
The origin of these deficiencies is now more or less understood from the investigations of 
the exact properties of the KS exchange-correlation potential V^ c s which the LDA or GGA 
does not satisfy: (i) V£ s (r) — > — 1/r as r — > oo for finite systems, PH|JTT[| (ii) the eigenvalue 



of the highest occupied single particle states is the negative of the ionization energy, [10112 



(iii) Vj£ s (r) is self-interaction free, and (iv) Vj£ s (r) exhibits an integer discontinuity upon 
the addition of an electron. [13,14] However, it has been recognized that approximations to 
the local KS potential obtained from the orbital dependent exact exchange energy functional 
through the so-called optimized effective potential (OEP) method [|1^,|1^] closely observe 
these conditions. [p~7| ,|T8|] 

In spite of these good features, the implementation of the OEP method is very com- 



plicated and its computational cost is excessively high. Only spherical atoms, flqJ17|J19 
jellium surfaces, [^tj jellium spheres, |HJ and periodic solids with a small unit cell [^] have 
been extensively studied, and there exists no investigation of molecular systems with the 
OEP method. However, about a decade ago, Krieger, Li, and Iafrate (KLI) proposed an 
approximation scheme to the OEP, and have shown that the KLI potential with the exact 
exchange-only energy functional (KLI-x) satisfies the above mentioned properties (i)-(iv) 
of Vj^ s , and accurately reproduces the OEP-x results for spherical atoms. |18| Within the 
KLI-x approximation, the workload in constructing the potential is significantly reduced, 
but this method is still relatively complicated compared with the LDA or GGA and studies 



are therefore limited. Atomic systems have been originally investigated by KLI, [M^ and 



more recently by Grabo et al. [33|] Bylander and Kleinman studied semiconductors with the 
KLI potential. |24]] Although molecules are a very good testing ground of any functional 
approximation, only some diatomic molecules have been investigated by Grabo et al. [23 

In this article, we report a study of several closed shell diatomic and polyatomic molecules 
within the KLI scheme. The questions we address are the performance of the KLI approxi- 
mation to the OEP with the exact exchange-only functional (KLI-x) and supplemented by 
LDA and Perdew-Burke-Ernzerhof GGA [25| correlation energy functionals (KLI-x+LDA-c 



2 



and KLI-x+GGA-c) for both structural properties and eigenvalue spectra. The comparison 
with and qualitative differences from the LDA/GGA and the Hartree-Fock (HF) method 
will be emphasized. In addition, since we incorporate the exact exchange, we can examine 
the accuracy of the LDA and GGA correlation energy functionals. The underlying com- 
putational algorithm is the three-dimensional (3D) finite-difference pseudopotential method 
|26| , |27|| and, to our knowledge, this is the first 3D KLI investigation of molecular systems. 

The organization of the paper is as follows. In Sec. O, the formulation of KLI is reviewed 
and the differences of the KLI-x from the HF method will be emphasized. The implementa- 
tion of KLI-x into our molecular code and its test results are described. In Sec. [TTT| , results 
and discussions for several molecules will be presented for both structural properties and 
eigenvalue spectra. Especially, comparisons between the KLI-x and the LDA/GGA or HF 
method are made, and physical origins of these differences will be discussed. In Sec. [TV], we 
conclude the paper by summarizing our work and future prospects. 



II. FORMALISM 

A. Approximation of OEP by Krieger, Li, and Iafrate 

We start from the spin-dependent KS equations for KS orbitals VVi( r ) for spin channel 
a =1,1 and KS eigenvalues e ai , with the possibility of fractional occupation numbers f ai 
(Hartree atomic units are used throughout the paper) 



i V 2 + ^(r) 



^o-i(r) = e CTi VVi(r), (1) 



with the spin density 



n rT {Y) 



E^IWr)! 2 , (2) 



1=1 



where N a is the number of occupied a spin orbitals. The effective KS potential (r) is 
composed of external, Hartree, and exchange-correlation contributions, 

V? S (v) = V ext (r) + V h (t) + V™{v), (3) 

where the exchange-correlation potential V^f(r) is given by 

KS = 5E xc [{n hni }] 

If the exchange-correlation functional is explicitly orbital-dependent £ , j BC [{n-[-,nj_}] = 
E X c[{i(j<Ti}}, the corresponding V^ff(r) = V^f p (r) can be obtained by the OEP integral 
equation [ff5l[16| 



E/<* / ^[^rCO-^COlCCOGoiCr'.r^Cr) + c.c. = 0, (5) 

ai 

where the quantity u xca i(r) is given by 
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MrJ_ Wr) • W 

The KLI approximation to the OEP is obtained by ignoring certain terms whose average 
over the ground-state spin density vanishes. Jl8| It can be also heuristically derived by 
approximating (e a j — e^) -1 , appearing in the Green's function in Eq. (|5|) by a constant 



j-j-l £<jj ^cji '-^Ca- 
in which case Eq. (|7]) reduces to integrals over the wave functions independent of the 
spectrum. Slater's approximation [p| can be obtained by simplifying the numerator of 
Eq. (0) further to i ) aj(' r ')'4 ) aj( r ) ~ ~~ r ')> thus ignoring — ^(r^^^r) term. |p~5|JT8 



Although the KLI approximation is generally applicable to any orbital dependent exchange- 
correlation energy functional, we will concentrate on the case of the exact exchange energy 
functional defined as 

ExiWai}] =--LL fnfaj j dv j dr | _ r /| ' ( 8 ) 

which has been actually used in this work. Here note that {VVi} are KS orbitals rather 
than HF orbitals. In this case, u XyCr i(r) is just the HF local exchange potential expression 
evaluated with KS (instead of HF) orbitals {ipai} 

u x ,ai(r) = — — — - J2 f^j( r ) K aji(r), (9) 

where 

K aji {r) = J dv J |r _ r/| . (10) 

The approximation of Eq. ([/]) leads to the KLI-x potential 

C L/ -*M = Vl(r) + £ ^U(V™-» - fl, iBi ), (11) 



1 n c 
i=i c 



where 



i= i n ^ 



is the Slater exchange potential, p8| with n ai {r) = f a i\tpai{ r )\ 2 - The quantities V x K J J i I ~ x and 
u X)C7 i are the expectation values of V^ li ~ x {y) and u x>ai {?) with respect to the KLI-x orbital 

ai, 

y^ti~ X = (lpai\V x KLI - X \Ai); U x>oi = (^K«n|V><n>- (13) 



The condition 



v£$r = ***** (14) 



is automatically satisfied with the exact exchange energy functional, |23| so only the remain- 
ing N a — 1 constants (V X K J^ I ~ X — u Xj0 -i) in Eq. ([Tlf) need to be determined by the solution of 
the following linear equation: 



N a -1 



X (<\n ~ M aji ){V™-* - u xM ) = (V x % - ,,,,„,). (15) 



i=l 

where 



n aj (r)n n 



M aji = dr U3K > alK \ i,j = l,...,N a -l. (16) 
J n a [r) 

The Slater potential V^ a (r) in Eq. flTT| ) ensures the correct asymptotic behavior of —1/r 



for large r when the exact exchange energy expression is employed. [16,18| It is known 



that V^f CT (r) alone is too deep, fTj| and is a factor of 3/2 larger than V^f = V^ D in 



the homogeneous-electron-gas limit. The second term in Eq. (fLip , which was originally 
missing in the Slater's proposition, enables Vj c J' I ~ x (r) to correctly reduce to V^ D in the 
homogeneous-electron-gas limit, and to preserve the property of the integer discontinuity. 
fl8l Also, due to Eq. (|T^) , the KLI-x highest occupied eigenvalues satisfy the following 
condition 

4n~ X = {^nM—V 2 + Vext + V H + V x KLI ~*)\^ Na ) 

= <VVivJ(-^V 2 + V ext + V H + u x , aNa )\iP aNa ) = eX (17) 

where e^' is the HF orbital energy expression evaluated with the KS (rather than HF) 
orbital Vw CT - 

At this point, we would like to emphasize that the KLI-x KS scheme yields a local 
potential, while the HF method gives a nonlocal exchange potential. It is well known that 
the HF method gives eigenfunctions and eigenvalues for the unoccupied orbitals which are 
inappropriate for describing actual excited states, because HF unoccupied orbitals do not 
have self-interaction corrections. This point will be further emphasized and discussed in 



Sec. II IB with specific examples. 



B. Computational Methods 

All the calculations in this work are based on the finite-difference pseudopotential scheme 
which we have previously applied to the study of neutral and charged water clusters. |27j 



In this formulation, quantities such as KS orbitals, densities, and potentials are expressed 
on grids in a rectangular simulation box, and the Laplacians or gradients are evaluated 
by higher-order finite-difference expressions. The most computationally demanding step 
in the construction of the KLI-x potential is the computation of the (N a — l)N a /2 pairs 
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/ dr / ^CTj( r OVVi( r O/l r — r '| i n Eq.@. Since these terms have the same structure as that of 
the Hartree potential with density n a {r') replaced by i/)* Ar')ip a i{r'), we use the techniques we 
have developed for the solution of the Poisson equation, and solve the discretized Poisson-like 
equations 

V 2 K aji {r) = -47r< i (r)V CTi (r). (18) 

To do so, we first generated boundary values for K a ji{v) at six sides of the simulation box 
using a multipole expansion. Multipole expansions have been performed until the boundary 
value error \5K a j i {r boundary )\/\K rT: j i {r houndar y)\ was less than 1CT 3 . Next, a coarse solution 
for K a ji(r) on the entire simulation grid has been generated by the fast Fourier trans- 
form method using a lower-order (3 points along each direction) finite-difference expres- 
sion. This solution has been subsequently relaxed by the iterative preconditioned conju- 
gate gradient method using a higher-order (typically 13 points along each direction) finite- 
difference method. Conjugate gradient relaxation steps have been performed until the error 
\8K a ji(r)\/\K a ji(r)\ is decreased to below 10~ 3 . For the preconditioning, a smoothing type 
preconditioner has been used. We observed that potential mixing rather than density mixing 
is more appropriate, which should be understandable considering that the KLI-x potential 
is orbital-dependent. 

In addition, for the GGA calculations, we have constructed GGA potentials with the 



scheme of White and Bird. [ 30 1 For other computational details including the implementation 



of the GGA potential, we refer to our other publications. |27|j3l|j32 



C. Test calculations 

To assess the accuracy of our 3D KLI-x implementation, we considered atoms and com- 
pared LDA and KLI-x results with those from an accurate ID radial atomic code. One test 
was for an Appelbaum-Hamann type local pseudopotential |33] modified to bind 5 electrons 
V(r) = -f erf(yfar) + {v x + v 2 r 2 )e~ ar \ where a = 0.6102, v x = 3.042, and v 2 = -1.732. 
This choice was made to reconcile the limitations of the 3D code that it can deal with 
only smooth pseudopotentials and the ID code that it can only produce accurate results for 
spherical atoms because of the central field approximation. Using a grid spacing h — 0.4 
a.u., simulation cell length L = 32 a.u., and finite-difference order 13, we obtained agreement 
with the radial code of < 10~ 3 a.u. /atom for both LDA and KLI-x total energies. For the 
hydrogen atom, we have used a local potential derived by Giannozzi |34| and again obtained 



an agreement of < 10~ 3 a.u. /atom using h = 0.4 a.u. For several other atoms considered 
in the next section (C, N, and O), we generated LDA pseudopotentials with the method of 
Troullier and Martins |$5|] and compared LDA total energies and eigenvalues from the two 
codes. With h = 0.25 a.u., the 3D code reproduced the ID atomic total energy with an 
accuracy of < 5 x 10~ 3 a.u. /atom which is sufficient for our purpose. 

Before closing this section, we would like to comment on the issues involved with the use 
of pseudopotentials. Although the rigorous and consistent procedure would be to employ 
pseudopotentials generated within the same functional approximations, we used LDA gener- 
ated pseudopotentials to carry out all other functional calculations for two reasons. First, we 
believe this procedure will not change the qualitative picture, since structural and electronic 
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properties are rather insensitive on the nature of the exchange-correlation energy functional 
that is used for the small Is cores of the first-row atoms considered here. Second, we 
note ambiguities associated with pseudopotential generation procedures. For example, we 
observed a small but sharp peak at the atomic center in the PBE GGA pseudopoential, 
which caused a serious convergence problem. |32| Moreover, in the KLI-x pseudopotential 
generation, there is an intrinsic problem of slow decaying tails due to the nonlocal nature of 
the exact exchange functional that somewhat influences the results. [0,0 These anoma- 
lous behaviors can be only alleviated by some post-processings which introduce additional 
arbitrariness in pseudopotenials. PH , |52"| 



III. RESULTS AND DISCUSSION 

A. Structural properties 

We first examine the structural properties of H 2 , CO, N 2 , H 2 0, and CH4 given in 
Table | and Table |J The binding energies Eg, equilibrium bond lengths r e , and vibrational 
frequencies u e of the diatomic molecules were determined from the total energy versus bond- 
length curves by fitting to the five parameter (D e , r e , uj e , a, and b) Hulburt-Hirschfelder 
function, |5S|J5|] 

V{r) = D e [(1 - e~ Px ) 2 + b(3 3 x 3 e- 2(3x {l + a(3x) - 1], (19) 

where x = r — r e and r e and D e are the equilibrium distance and bonding energy, using 
the simulated annealing method. jHJ The geometries of the polyatomic molecules have been 
optimized by employing Hellmann-Feynman quantum forces. A grid spacing h = 0.25 a.u. 
has been used for all the molecular calculations. For the calculation of binding energies, the 
energies of each pseudoatom and molecules have been calculated by the same method and 
with the same h. We expect this procedure will result in a systematic cancellation of errors, 
and in addition we may treat the atom without central field approximation. 

Our LDA and GGA results are essentially a reproduction of previous studies. P^H ^l 



The LDA consistently overestimates the stability of the molecules and the GGA typically 
cures this tendency of overestimation to a large extent. The HF approximation on the 
other hand substantially underestimates the binding energy. p9|M4lj Bond lengths calculated 
by the LDA are usually overestimated and one obtains corresponding underestimations of 
vibrational frequencies for the diatomic molecules. The HF method gives the opposite 
behavior. The KLI-x approximation typically gives similar results as HF for the structural 
properties which are mostly determined by the total energy functional. 03] It is well known 



that the HF method gives the wrong sign of the dipole moment of CO, and the LDA corrects 
this. Again, the KLI-x calculation yields the wrong sign of the dipole moment of CO in 
agreement with the HF method. The difference between magnitudes of the KLI-x and HF 
dipole moment of CO results from the fact that the HF value has been evaluated at the 
experimental bond length. Our KLI-x calculation at the experimental bond length gave 
-0.275 Debye. 

The addition of an LDA or GGA correlation energy functional to the exact exchange 
energy gives mixed results: Binding energies are increased, hence improved the KLI-x values, 
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but bond lengths are decreased and become worse. So, we can conclude that the two 
correlation functional approximations do not give systematic improvements over the exact 
exchange energy functional, which indicates that these approximations are not accurate 
descriptions of correlation effects in molecules. 



B. KS Eigenvalues and Excitations 

In Fig. |l|, we compare the LDA-xc and KLI-x potentials of H 2 to reemphasize the 
different nature of the exchange (-correlation) potentials from the two methods. While the 
LDA-xc potential decays exponentially, the KLI-x potential reproduces the correct —1/r 
asymptotic behavior, and this causes marked differences in the eigenvalue spectra. This 
—1/r decay property lets the high-lying KS unoccupied eigenvalues correctly approach a 
Rydberg series, and is also important for other properties that are sensitive to the outer 
part of the charge density, such as polarizabilities. 

Eigenvalues of the highest occupied orbitals of the molecules considered in the previous 
section at their equilibrium geometries (Tables [I] and [TT|) are listed in Table |TJ. Compared 
with the LDA or GGA, the KLI-x values agree very well with experimental ionization po- 
tentials, and they are very close to HF values. It is well known that HF gives a good 
approximation to the ionization potential via Koopman's theorem due to the fact that omis- 
sion of correlation tends to be cancelled by the neglect of relaxation in the "frozen orbital" 
approximation. Q| On the other hand, in the exact DFT, the highest occupied orbital 



eigenvalue equals the ionization potential without any relaxation correction [[Kjillil as we 
typically see in our KLI-x calculations, hence, in principle DFT should yield highest occu- 
pied orbital eigenvalues in better agreement with experimental ionization potentials than 
HF. However, in the usual LDA or GGA, the highest occupied eigenvalues are in serious 
errors due to the wrong decay property of exchange-correlation potential exemplified in Fig. 
[l]. The errors are about 0.2 a.u. for the molecules we have studied as shown in Table [HI. 



On the other hand, due to Eq. (|TT|) , KLI-x should give highest occupied orbital eigenvalues 
in good agreement with experimental ionization potentials when HF values also agree with 
experimental ionization potential. The addition of LDA or GGA correlation potentials to 
the KLI-x potential usually lowers KLI-x highest occupied orbital eigenvalues, but this leads 
to too large ionization potentials. 

Among the molecules examined in this work, N 2 is a particularly interesting case which 
shows the advantage and disadvantage of the KLI-x scheme at the same time. First, note 
that the HF results are in qualitative disagreement with experiment: it puts the lir u state 
higher than 3a g state as the highest occupied orbital in contrast to the experimental data. 
This incorrect ordering of the first two ionization potentials is a well-known example of the 
breakdown of the HF picture and indicates a large effect of correlation upon the orbitals in 
N 2 . MM On the other hand note that KLI-x gives the correct ordering of the two orbitals, 



which implies that the local nature of the KLI-x potential can make a difference even for the 
occupied-orbital eigenvalue spectrum. The fact that the LDA also gives the correct ordering 
(although the value is worse than the KLI-x one due to the above-mentioned reasons) shows 
that there are cases where the local potential in the KS scheme is superior to the nonlocal 
HF potential. The highest occupied state in the KLI-x scheme e^ g I ~ x is equal to e^f 
[according to Eq. (|T7|)] which results in the biggest KLI-x error among the molecules studied 



S 



since this also has the largest HF error. Better quantitative agreement with experiment can 
be obtained by incorporating the correlations, e.g. through many-body perturbation type 
approach. |44] 



In addition to the eigenvalues of the highest occupied orbitals, we have also examined 
eigenvalue differences between the highest occupied and the lowest unoccupied orbitals, 
which are compared with experimental lowest triplet and singlet vertical excitation energies 
in Table |TV|. First, we point out that the HF calculations give positive orbital energies for 



all the virtual orbitals in the molecules considered here. |39| , |44|1 This is because unoccupied 
orbitals in HF do not really correspond to excited states of the system, in which the excited 
state electron would have been removed from one of the lower states and acted on by N — 1 



remaining electrons, |29] but rather the bound states (if any) of the singly- charged negative 
ion, in which the extra electron sees the field due to other N electrons. [B9| To make this 
argument more explicit, we follow Slater [^5| and rewrite Eqs. (|9|) and (|10"D (now using the 



HF orbitals ip^f instead of KS orbitals i[) ai , so u^i is the true HF local exchange potential) 



as 



HF( , f , , [gfe / CTJ ^(r)< F *(rQ^(r)^(rQ]/[^(r)^(r)] 
u x , a A r ) = ~ J dr ' ( °) 

and identify 

Ef=i^^(rX^(rQ^(r)^(rQ 

^ F *(r)€f(r) [Zi) 

as an exchange charge density. Because the ip^ F, s are orthonormal, the exchange charge 
density integrated over dr' will be minus one for the occupied orbital (j = i term exists in 
the summation J2j=i an d f a % = 1), but zero for the unoccupied orbital (ip^ is not included 
in J2j=x)- Hence, the self-consistent field in which an unoccupied orbital moves is that of 
the nucleus and all N electrons instead of N — 1 electrons. 

On the other hand, in DFT, as mentioned in Sec. |, differences of KS eigenvalues are well- 
defined approximations for excitation energies. [0] Umrigar et al. have confirmed that this 
is the case for several atomic systems by showing that the exact KS eigenvalue differences 
obtained from KS potentials derived from accurate QMC densities almost always lie between 
the triplet and singlet experimental excitation energies. || They further claimed that this 
is because of the fact that, in addition to the well known exact properties of V^ c s (Sec. |), 
Vj^ s agrees with the quasiparticle amplitude |L0| for r — > oo not only up to 1/r but up to 



order 1/r 4 inclusive. However, in the LDA or GGA, no bound unoccupied state exists in 
many cases, due to the shallow nature of the corresponding exchange-correlation potentials. 
We can expect that the KLI-x is a much better approximation than the LDA/GGA, since 
it satisfies Eq. (EHl) and has the correct asymptotic —1/r behavior. Our KLI-x results 
indeed always give negative virtual orbital eigenvalues for the studied molecules and also 
show good agreement with experimental values. The single exception is H 2 in which the 
difference between the highest occupied and lowest unoccupied orbital eigenvalues is bigger 
than the experimental triplet and singlet excitation energies by ~ 0.07 a.u. Further detailed 
agreements with different multiplet states can be obtained by employing a more involved 
theory such as DFT perturbation theory W\M or time- dependent DFT. HJBI] 
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We close this section by reconsidering the meaning of eigenvalues in the OEP-x (or KLI- 
x) and HF schemes to make the differences between the two be clear. In the fractional 
occupation number formalism, KLI have shown that 

dE 1 

Of. = (^l(-2 V2 + v ** + V H + = t H J\ (22) 

where e^ p> is the HF orbital energy expression for the orbital ai evaluated using the OEP-x 



orbitals ip a {. []T8j We note that this can be generalized to the correlated OEP case. |18|j46 
Equation (|22| ) can be rearranged to the form in which the relationship between the two 
methods is transparent: 

£ EP ~ X + (VVil ( u x,ai ~ V x OEp -n\^) = e»J'. (23) 

Note that from the condition V X E ^~ X = u x ^N a , satisfied by both OEP-x and KLI-x [Eq. 
(PDl, Eq. fllTD is reproduced from Eq. ( p3| ) for the aN a state. Hence, the eigenvalue 
difference between the highest occupied orbital state ip a N a and the unoccupied orbital state 
ip aa is given by 

= (24) 

Assuming that the differences between the HF and OEP-x orbitals are negligible, we can 
see that HF orbital eigenvalue differences are much bigger than OEP-x values, because the 
(i J aa\u x ,a-a\' l Paa) term is small due to the property of the exchange charge density [Eq. (pip ] 
for the unoccupied orbital as discussed above, while the (ip a a\V x 0EP ~ x \i J aa) term is clearly 
quite large. 



IV. CONCLUSIONS 

In this work, we have studied structural and excitation properties of small molecules 
using the KLI approximation to the OEP with the exact exchange energy functional only 
(KLI-x) and augmented by LDA and GGA correlation functionals. For structural properties, 
the KLI-x gave comparable results to those of the HF method. For excitation properties, 
the KLI-x results in good eigenvalue spectra for both the highest occupied and unoccupied 
orbitals, since it has the correct asymptotic large r behavior of Vj£ unlike the LDA or 
GGA. Especially, we find that unoccupied orbital eigenvalues are better described by KLI-x 
than by HF, which illustrates the advantage of the KS scheme in general. 
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TABLES 



TABLE I. Binding energies, bond lengths, and vibrational frequencies of H 2 , CO, and N 2 , 
calculated with different methods. In addition, the dipole moment of CO has been given, whose 
sign has been defined such that C~0 + is positive. Energies and bond lengths are given in a.u., 
vibrational frequencies in cm -1 , and dipole moments in Debye. 







H 2 






CO 






N 2 






Eb 


r e 




E B 


r e 






Eb 


r e 




LDA-xc 


0.184 


1.44 


4228 


0.466 


2.13 


2193 


0.248 


0.396 


2.09 


2385 


GGA-xc 


0.168 


1.41 


4228 


0.424 


2.12 


2166 


0.265 


0.368 


2.09 


2383 


KLI-x 


0.136 


1.39 


4647 


0.287 


2.06 


2445 


-0.163 


0.176 


2.01 


2634 


+LDA-C 


0.185 


1.36 


4765 


0.342 


2.05 


2595 


-0.184 


0.271 


2.00 


2710 


+GGA-c 


0.166 


1.38 


4734 


0.361 


2.06 


2641 


-0.170 


0.273 


2.00 


2805 


HF-x a 


0.134 


1.39 


4582 


0.293 


2.08 


2431 


-0.279 b 


0.195 


2.01 


2730 


Exp. a 


0.174 


1.40 


4400 


0.414 


2.13 


2170 


0.123 


0.364 


2.07 


2358 



a Reference |3l| . 

b At the experimental bond length. 



TABLE II. Binding energies and bond lengths of CH4 and H 2 0, and the bond angle of H 2 0. 
Binding energies and bond lengths are given in a.u., bond angles are given in degrees. 







CH 4 




H 2 






E b 


r e (C - H) 


E b 


r e (O - H) 


9 (H-O-H) 


LDA-xc 


0.743 


2.07 


0.420 


1.82 


104.2 


GGA-xc 


0.681 


2.06 


0.370 


1.81 


104.5 


KLI-x 


0.536 


2.05 


0.257 


1.76 


106.0 


+LDA-C 


0.678 


2.04 


0.344 


1.76 


106.0 


+GGA-c 


0.672 


2.03 


0.339 


1.77 


105.3 


HF-x a 


0.522 


2.05 


0.258 


1.78 


106.1 


Exp. a 


0.668 


2.05 


0.374 


1.81 


104.5 



^Reference |3Sj. 



13 



TABLE III. Absolute values of highest occupied orbital eigenvalues of the H%, CO, N2, H2O, 
and Ci/4 at their equilibrium geometries (Table I and II). For N2, absolute values of the next 
highest occupied orbital eigenvalues are also given. Note the differet ordering of HF and KLI-x for 
this case. Experimental values are ionization potentials. Energies are in a.u. 





H2 (-eia 9 ) 


CO (-e 5a ) 


^2 (-e 3 <x g / - eivrj 


H 2 (-e 16l ) 


CH4 (— eit 2 ) 


LDA-xc 


0.37 


0.34 


0.38/0.44 


0.26 


0.34 


GGA-xc 


0.38 


0.34 


0.38/0.43 


0.27 


0.34 


KLI-x 


0.60 


0.55 


0.64/0.69 


0.50 


0.54 


KLI-x + LDA-c 


0.60 


0.61 


0.69/0.75 


0.56 


0.60 


KLI-x + GGA-c 


0.58 


0.59 


0.67/0.74 


0.55 


0.58 


HF-x a 


0.60 


0.55 


0.64/0.62 


0.51 


0.55 


Exp. a 


0.58 


0.58 


0.57/0.62 


0.46 


0.53 



a Reference |3El. 



TABLE IV. Differences between highest occupied and lowest unoccupied orbital states eigen- 
values of H2, CO, N2, H2O, and CH4 at their respective equilibrium geometries given in Table I 
and II. A — mark indicates that the lowest unoccupied orbital state is in the continuum and hence 
unstable. Experimental values are vertical excitation energies to final triplet and singlet states 
shown in parentheses. Energies are in a.u. 





H2 (eio-,, - eio-J 


CO (e 27r - e 5(T ) 


^2 (ein g -e3a 9 ) 


H2O [€4 ai — €lf )1/ 


CH4 [€3 ai — e 


LDA-xc 




0.26 


0.30 






GGA-xc 




0.26 


0.31 






KLI-x 


0.47 


0.28 


0.34 


0.34 


0.42 


KLI-x + LDA-c 


0.50 


0.28 


0.34 


0.36 


0.46 


KLI-x + GGA-c 


0.49 


0.27 


0.34 


0.37 


0.45 


HF-x 












Exp. a 


0.42 ( 3 S+) 


0.23 ( 3 n) 


0.28 ( 3 S+) 


0.26 ( 3 Si) 


0.40 ( 3 T 2 ) 




0.46 CU U ) 


0.31 (^n) 


0.34 Cn g ) 


0.28 ^Bi) 


0.41 (!r 2 ) 



1 Reference [3S]. 
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FIGURES 

FIG. 1. Comparison of the LDA exchange-correlation potential and the KLI exchange potential 
for H2 with H atoms located at ±0.7 a.u. Note that the LDA exchange-correlation potential 
incorrectly decays exponentially, while the KLI exchange potential decays with — 1/r as in the 
exact Kohn-Sham exchange-correlation potential. 
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